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Abstract. In this paper we study the large time behavior of solutions to a 
nematic liquid crystals system in the whole space R 3 . The fluid under consid- 
eration has constant density and small initial data. 



1. Introduction 

In this paper we consider the asymptotic behavior of solutions to the simplified 
model of nematic liquid crystals (LCD) with constant density: 

u t + u ■ Vu + Vp = vAu - V • (Vd <g> Vd), 
(1.1) d t + u-Vd = Ad- f(d), 

V -u = 0. 

The equations are considered in R 3 x (0,T). Here p : R 3 x [0,T] ->• K is the fluid 
pressure, u : R 3 x [0, T] -> R 3 is the fluid velocity and d : R 3 x [0, T] ->■ R 3 is 
the direction field representing the alignment of the molecules. The constant v > 
stands for the viscosity coefficient. Without loss of generality, by scaling, we can set 
v = \. The force term V<i ® Vd in the equation of the conservation of momentum 
denotes the 3x3 matrix whose ij-th entry is given by "V^d • Vjrf ' for 1 < i,j < 3. 
This force Vd ® Vd is the stress tensor of the energy about the direction field d, 
where the energy is given by: 

1 

2 

where 



/ |Vd| 2 d:r+ / F(d)da; 
Jk 3 Jr 3 



M 3 



F(d) = — (\df - l]\ /(d) - VF{d) = -(|d| 2 - l)d, 

for a constant r] in this paper. We note that F(d) is the penalty term of the 
Ginzburg-Landau approximation of the original free energy of the direction field 
with unit length. 

In this paper we consider the following initial conditions: 

(1.2) u(x,0) = u (x), V-u o = 0, 

(1.3) d{x,0) = d (x), |db(x)|=l, 
and 

(1.4) uoGi^OR 3 ), do - w Q e H 2 (R 3 ), 
with a fixed vector wq € S" 2 , i.e., liOol = 1- 
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The flow of nematic liquid crystals can be treated as slow moving particles where 
the fluid velocity and the alignment of the particles influence each other. The 
hydrodynamic theory of liquid crystals was established by Ericksen [7] [5] and Leslie 
[TBI [17] in the 1960's. As Leslie points out in his 1968 paper: "liquid crystals are 
states of matter which are capable of flow, and in which the molecular arrangements 
give rise to a preferred direction" . There is a vast literature on the hydrodynamic 
of liquid crystal systems. For background we list a few, with no intention to be 
complete: [illlEllIllinillHlIIlllllllllillllE]- In particular, the asymptotic 
behavior of solutions to the flow of nematic liquid crystals was studied for bounded 
domains in )19i 28 j. It was shown in [28] that, with suitable initial conditions, the 
velocity converges to zero and the direction field converges to the steady solution 
to the following equation 



In [28], Lemma 2.1 the Lojasiewicz-Simon inequality is used to derive the con- 

vergence when fl is a bounded domain. Lack of compactness considerations do not 
allow us to use similar arguments in the whole space R 3 . 

In this paper we consider the asymptotic behavior of the solutions to (|l.ip in 
the whole space M 3 . subject to the additional condition on the direction field which 
insures that the direction tends to a constant unit vector wq, as the space variable 
tends to infinity: 



This simplifies the situation and allows us to obtain the stability without need- 
ing the Liapunov reduction and Lojasiewicz-Simon inequality, since wq is a non- 
degenerate steady solution to (|1.5[) . 

We start from the basic energy estimates (|2.16[) and Ladyzhenskaya estimates 
(|2.17|) [15] [6] (see the extension to the whole space in appendix of this paper) for 
the system (|l.ljl . We then establish the convergence of the direction field d to 
the constant steady solution wq based on Gagliardo-Nirenberg interpolation tech- 
niques. More precisely, the convergence obtained is in L P (M 3 ) for any p > 1, with 
an algebraic decay rate of (1 + t) - ^ 1- ?). Wc then focus on the conservation of 
momentum equation in (|1.1[) . We apply the Fourier splitting technique [2"51 [2^1 |2"T] 
to obtain the decay of the velocity u with an algebraic decay rate of (1 + t)~i in 
L 2 (R 3 ) norm. This rate coincides with the decay rate of solutions to Navier-Stokes 
equations with a force decaying at a rate of (1 + i) - * [23] . 

The existence of global regular solutions of ([l.ip with the initial and bounadry 
conditions has been established in |19] (in [B] for nonconstant density) provided 
that the viscosity constant is large enough or initial data are small enough. Based 
on the arguments in [BJ the existence of global regular solutions of (|1.1[) . for small 
initial data, is established in the appendix as follows: 

Theorem 1.1. Let uo and do satisfy U.2)) - (L4\ ). Assume that uq s H 1 (S. 3 ) and 
da — wq G i7 2 (R 3 ) n L 1 (R 3 ) for a unit vector wq. There is a positive small number 
eo such that if 



(1.5) 





(1.6) 



lim do(x) = wq. 

\x\—^oo 



(1.7) 



IKIIffi(K3) + \\<k - w \\ 2 H 2 {R3) < e 
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then the system 1 1.1]) has a classical solution (u,p,d) in the time period (0,T), for 
all T > 0. That is, for some a € (0, 1) 

u e C 1+a/2 ' 2+a ((0,T) x M 3 ) 

(1.8) VpeC a/2 ' a ((0,T)xR 3 ) 

d g C* 1+q/2 ' 2+q ((0,T) x R 3 ). 

And the solution (u,p,d) satisfies the following basic energy estimate and higher 
order energy estimate (also called Ladyzhenskaya energy estimate in [B] and |19j ) 

(1.9) / M 2 + |Vd| 2 + 2F(d)dx + 2 f [ |Vu| 2 + \Ad - f(d)\ 2 dxdt 
Jr. 3 Jo Jr 3 

< ll w o|li2 (R3) + ||Vdo|li2( R 3) 

(1.10) / |Vu| 2 + \Ad\ 2 dx+ [ [ \Au\ 2 + \\7Ad\ 2 dxdt 
Jr 3 Jq Jr 3 

< Cdluoll^i^) + \\d - wq\\ 2 H 2 (r3) ). 
Furthermore, the solution d satisfies 

(1.11) / \d{t) - w \dx < {C t + | d - w \dx)e ct 
Jr 3 Jr 3 

with the constants Co and C depending only on initial data and on n, respectively. 
For the smooth solution obtained in Theorem 11.11 our main asymptotic result 

is: 

Theorem 1.2. Let (u,p,d) be smooth solution obtained in Theorem \l.l\ Assume 
additionally uq G L 1 (IR 3 ) and do — wq G L p (M. 3 ), for any p > 1 and a unit vector 
wo ■ There exists a small number eo > such that if 

(!- 12 ) IKHfrifRS) + ||d - woWh'QP) < e o. 

then 

(1.13) \\d(-,t) - w \\ LPm < C\\d - ttf ||iP(R»)(l + t)-^ 1 -^, 

(1.14) \\V(d(-,t)-wo)\\h m <C(l + t)-i, 

(1.15) \n;t)\\h { R 3 ) <C(l + t)-i, 

where the various constants C only depend on initial data. 

The paper is organized as follows: in Section 2 we establish the decay for the 
difference d — wo , using the basic energy estimate (|1.9[) and the Ladyzhenskaya en- 
ergy estimate (|1.10|) . Combining the decay of d — wo and Fourier splitting technique 
[23] . in Section 3 we obtain an algebraic decay for the velocity u in L 2 (E 3 ). In the 
appendix, we sketch a proof for the existence theorem ll.il 



4 



MIMI DAI, JIE QING, AND MARIA SCHONBEK 



2. Convergence of the direction field 

In this section we study the L p decay of the direction field d — wo and the decay 
for the first derivative. The first step is to derive a uniform estimate of d — wq in 
L 2 (R 3 ). This yields a uniform estimate for d — wo in L P (R 3 ) for any p > 1. This L p 
estimate, is the basis to establish the decay results. For ease of reading we state the 
basic energy estimate and Ladyzhenskaya energy estimate satisfied by the smooth 
solution [u,p,d) (see appendix for details), 

rT 

f ||Vrf||| 2 + \\F< 
<|| Mo ||| 2 + ||Vd || 2 



(2.16) \\u\\U + ||Vd||£ a + \\F(d)\\ L i + 2 / ||V«||£ a + ||Ad - f(d)\\i 2 dt 
\h + \\Vd \\l 2 



(2.17) ||V U || 2 i2 + ||Ad||| 2 + f T \\Au\\l 2 + \\VAd\\l 2 dt <C(\\u \\l 1 + \\d -w \\ 2 H2 ). 

Jo 

In the sequel we need to use a Gagliardo-Nirenberg interpolation inequality. For 
completeness we recall from QT] the inequality here 

Proposition 2.1. [IT] Let w € W m ' p (M n )nL''(M. n ), for 1 < p < oo and 1 < q < oo. 

Then 

(2-18) \\D k w\\ L r (Mn) < C\\D m w\\l P(Rn) \\w\\l~^ 

for any integer k £ [0, m — 1], where 

1 k ,1 m . , ,1 

2.19 - = - + a( + (l-o)- 

r n p n g 

wif/i a 6 [~, 1], eii/ier if p= \ or p> \ and m — k—- J\fU {0}, wMe a G 1), 
i/p > 1 and m - A: - - £ TV U {0}. 



2.1. Uniform estimate of d — wq in i 2 (R 3 ). In this part, we show that the in- 
tegrals 

d(x, t) — u>o\ 2 dx and / / \\7(d(x, t) — wo)\ 2 dxdt 

JO Jr 3 

are uniformly bounded by the initial data, applying estimates for the solutions 
obtained in Theorem ll.il We have the following lemma, 

Lemma 2.2. Let d be the solution obtained in Theorem ] 1. 11 There exists eg suffi- 
ciently small so that if || wo || z, 2 (M 3 ) + II^ 7 '^o||l 2 (r 3 ) < e o- Then 

|2j„ t / / ±\ „,,_M2, 



(2.20) / \d(x,t)-w \ 2 dx+ / / |V(d(x,i) - w )\ 2 dxdt < C, 
Jr 3 Jo Jr 3 

where the constant C depends on the initial data and the norm \\d$ — wo||i2( R 3). 

Proof: Since |wo| = 1 and f(wo) = ^{\wo\ 2 — l)wo = 0, the second equation in 
(jl.lj) can be expressed as 

(2.21) (d - w )t + u-\7{d- w ) = A{d - w ) - f(d) + f(w ). 
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Applying the mean value theorem for vector valued functions, we have 

(2.22) f(d) - f(w ) = (J Df(w + s(d - w ))d^j ■ (d - wo), 

where Df denotes the Jacobian matrix of /. Multiplying (|2.21[) by d — wo yields 

(2.23) ~ f \d-w \ 2 dx 

z at j R3 

[u ■ V(d — Wo)](d — wo)dx 

+ I A(d — wo) ■ (d — wo)dx 

i 

(d - wo) T I / Df(w + s(d - w ))ds I (d - w )dx, 
! Vo / 

= -h + h + h- 
The three terms I4, I§ and 1$ are estimated as follows: 

(2.24) \h\ = \ lim / [u ■ V{d - w )](d - w )dx\ 

— I \d — wo\ 2 u ■ nda\ 

2 JdB R 



= I lim - 

_R^oo 2 



<C lim ( f \d-wo\ 4 da) ' ( [ \u\ 2 da 

R ^°° \JdB R J \JdB R 

<C*lim (f \d-wo\do] ([ \u\ 2 da 

\JdB n J \JdB R 



Denote the set A = |i? : J 9Br \d — wo\da > m|, for a certain constant M > 0. A 
is closed and its complement A c is open. Write 

/ \d(t) -w \dx = [ f \d(t) - w \dadR. 

JR 3 JO JdB a 

Recall that by Theorem the last expression is bounded by the initial data, for 
any fixed time t > 0. Thus by Chebyshev's inequality we have 

(2-25) 

where \x denotes the measure of a set, and C denotes a constant depending only on 
the initial data. Since the energy estimate (|2.16[) implies that 



[ \u\ 2 dx= [ [ \u\ 2 dadR 

JR 3 JO JdBn 



is bounded by initial data, there exists a sequence {Ri}°^ 1 C A c with Ri — > 00 such 
that 

(2.26) / \u\ 2 da -> 

JdB R . 



Combining the inequalities ff2T24]) . (f2T25|) and ([236]) . yields for all £ > 0, 
(2.27) h = 0. 
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For .Z5, we have 

(2.28) 7 5 =lim[/ 9{ * d ~ Wo) ■ (d - w Q )dcr - I \V(d - w )\ 2 dx] 



R ^°° JdB R on j Br 
The boundary term is estimated as follows: 

(2.29) / ^^.(d- W0 )da 
JdB R on 

<([ \V(d-w )\ 2 da] 1 ([ |d-w | 2 do-V 

\JdB R ) \JdB R J 

<c(f \V(d-w )\ 2 da) ' ([ \d-w \da) ' 

\JdB R J \JdB R J 

From Theorem ll.il we have that || V(d — wo)||l 2 (m 3 ) is uniformly bounded for any 
fixed t > 0, \\d(t) — woIIl^r 3 ) has a time dependent bound. For the inequality 
(|2.29l) . we apply a similar argument used to derive (|2.27p and obtain the existence 
of a sequence Ri approaching infinity satisfying 

f d(d-w ) , u „ 

Inn / • (d — w )da — > 

R i^°°JdB R . on 

It follows then that 

(2.30) h = - [ \V(d-w )\ 2 dx. 



With respect to I§, let d = wq + s(d — wq) and z — d — wq. Using the definition 
of f(d) = ^(\d\ 2 — l)d. a straightforward calculation yields 

(d-w ) T (J Df{w + s(d-w ))ds\ (d-w ) = z T J Df(d)dsz 
= f\(2d 2 + \d\ 2 - \)z\ + (2d 2 + \d\ 2 - l)z 2 + (2d 2 + \d\ 2 - l)z 2 



o 



*3 

4did 2 z 1 z 2 + 4did 3 ziz 3 + Ad 2 d 3 z 2 z 3 ]ds 



= ( [2d 2 1 z 2 + 2dlz 2 + 2djzl + 4d 1 d 2 z 1 z 2 + 4:d 1 d 3 z 1 z 3 
Jo 

+ Ad 2 d 3 z 2 z 3 ] + (\d\ 2 - l)(z 2 + z 2 + zfjds 

= f 2[d 1 z 1 +d 2 z 2 + d 3 z 3 } 2 + (\d\ 2 ~l)\z\ 2 ds 
Jo 

= [ 2[d-z] 2 + (\d\ 2 -l)\z\ 2 ds. 
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In the above equation, the third and forth equality comes from regrouping terms 
and completing a perfect square. Thus, Iq can be written as 



1 



i 

2 



(2.31) h = — 2/ / 2[(w + s{d-w ))-(d-w )] 

V JM 3 JO 

+ (\wq + s(d — wo)\ 2 — l)\d — wo\ 2 dsdx 

= 7? [ 2[w ■ (d - u) )] 2 + 3w ■ (d - wo)\d - w \ 2 + \d — w \ A dx 

T Jr 3 

If 3 1 

= 5- / 2ku ■ (d - w ) + -\d - w Q \ 2 } 2 --\d- w \ A dx 

V Jr 3 4 8 

- 7T~2 / \d-w \ 4 dx. 
Combining (f2~2"31 and the inequalities P^7|) . (|2~3U1) . and (|2~BT1) gives 

(2.32) ~ [ \d~w \ 2 dx+ [ |V(d- w )\ 2 dx < / |d-u; | 4 da;. 



The right hand side of the inequality (|2.32[) can be estimated as 
\d-w a \ A dx<{( \d~w \ 2 dx) 1/2 ( [ \d~w \ 6 dx) 1/2 

<C([ \d-w \ 2 dx) 1/2 ( f \V(d-w Q )\ 2 dx) 3/2 

JR 3 JR 3 

<C [ \d-w Q \ 2 dx( [ |V(d- w )| 2 ^) +7^/ |V(d-u; )| 2 dr. 

JR3 Vr 3 / 2 J R 3 

Gagliardo-Nirenberg interpolation inequality (Proposition 12 . ll) yields 
||V(d - w )||i2 (R 3) < C||d - w ||L 2 (R3)||A(d - w )||l 2 (r3)- 
Combining the last two inequalities with (|2.32[) gives 

(2.33) — : / \d-w \ 2 dx+ [ \\7(d-w )\ 2 dx 

">t Jr3 J R 3 

<c([ \d~w \ 2 dx^) [ \A(d-w )\ 2 dx. 
\Jr 3 ) Jr 3 

Denote 4>(t) = J R3 \d(t) — w \ 2 dx. Then 

^ <C [ \A(d-w )\ 2 dxdt. 

<P Jr 3 

Integrating the last inequality over [0,t] yields 



1 1 



</>(t) 0(0) " Jo 
Thus, 

0(0) 



< C ( I \A(d~w )\ 2 dxdt. 
Jo Jr 3 



l-CmloI R 3 \Md-w )\ 2 dxdt 



8 MIMI DAI, JIE QING, AND MARIA SCHONBEK 

From the basic energy estimate (|2.16[) . and the hypothesis we have 

/ / \A(d-w )\ 2 dxdt< IKHia + HVdollia <eo- 
Jo Jr 3 

Assume that e is so small that Ceo0(O) < 1/2, then 

CV>(0) / / |A(d- w )\ 2 dxdt < 
Jo Jr 3 2 

Hence for any t > 0, 

m < 20(0), 

that is, 

(2.34) / \d{t) -w \ 2 dx < 2 / \d -w \ 2 dx. 

Jm. 3 Jm 3 

Due to the estimates (|2.34[) and (|2 .33[) . we have 



d 
It 



\d-w \ 2 dx+ \V(d~w )\ 2 dx <C I |A(d-io )| dx, 

JR 3 



where the constant C only depends on the initial data. Integrating over [0, t], by 
the basic energy inequality (|2.16[) it follows that 



\d(t) -w \ 2 dx+ [ [ \V(d-w )\ 2 dxdt 
Jo Js. 3 



< I \d Q -w \' z dx + C I I \A(d-w )\ 2 dxdt 

< C. 



Thus, 



(2.35) / / |V(rf— w )\ 2 dxdt < C, 

Jo Jr 3 



where the constant C only depends on initial data. This completes the proof of the 
lemma. 

□ 

The following auxiliary estimate shows that, provided the initial data is small 
enough, the norm ||cf(-,£) — wo\\l°°(r 3 ) wm be as small as necessary. This smallness 
yields that \d\ will be close to 1, for all time. 

Lemma 2.3. Let d be the solution obtained in Theorem ] 1. 11 Then 
(2.36) \\d(;t) -tltolU-CRa) < CUVdWII^a^A^i)!!^), 

where C is an absolute constant. 
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Proof: Let T denote the Fourier transform. By Lemma 12. 2\ we can take the 
Fourier transform of d — wq 

(2.37) 

\\d(;t)-W \\ L oo m < [ \T{d-w M 
JR 3 

= ! nd-woM+ [ \F(d-w M 

= / 4 • \Z\\F(d ~ woM + I l4 ■ \?WHd - wM 

J\e,\<\ m J\t\>\ \£r 

^([ m2 d ok[ mnd-w )fd& + ([ jLdoHf \end- Wo )\ 2 d& 

<C(/ ^dr)i([ \\7(d-w a )\ 2 dx)i +C([ ^dr)?( / | A(rf - Wo )| 2 dx) ^ 

Jr<A r JR 3 Jr>A r JR 3 

< CA*||Vd(i)|| L 2 (R3) + CA-'||Ad(i)|| L 2 (H3) . 

To find a A that minimizes the right hand side of the last inequality, take the 
derivative in A and set the right hand side equal to zero: 

\i\\Vd(t)\\ L * m = A-*||Ad(t)|| La(B s), 

yielding 

A = || Ad(t) ||i3( B 3)/|| Wd(t)\\ L 2(^3y 

Thus using this A in f)2 .37[) gives inequality (|2.36|) . and the proof of the Lemma is 
complete. 

□ 

Corollary 2.4. Suppose the initial data ||ito Hjt 1 + \\do — wo\\h 2 are small enough. 
Then, \d(x,t)\ > \. 

Proof: It follows combining ([235]) . ([2~P7|) . (|2~3^1) since \wq\ = 1. 

□ 



2.2. Uniform estimate of d — wo in L p with any p > 1. Here we show provided 
the data is small enough, all the L p norms of d — wq are bounded. 

Lemma 2.5. Let d be the solution obtained in Theorem ] There exist X p , de- 
pending on p, so that if \\uo\\%it M 3) + \\do — w o||jj2( H 3) < \, then for p > 1 

(2.38) -[ \d(x,t) - w \ p dx+ ^ / / |V|d(s,t) - w | p/2 | 2 da;di 

PJr3 V Jo Jr3 



< C p | do - u>o| p dE, 

JR 3 

where the constant C p depends on p and X p . And for p = 1 we ftawe 

/ |d-w |dr< / |d -wo|dx+ / / |V(d - w )| 2 + |A(d - w )| 2 cfe(ft 

JR 3 JR 3 JO JR 3 

< c. 
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Proof: Recall that since f(wo) = 0, we have that the direction equation can be 
rewritten as 

(d - w ) t + u • V(d - w ) = A(d - w Q ) - /(d) + f(w ). 

Multiplying the last equation by (d — Wo)\d — w \ p ~ 2 , for any p > 2 (or alternatively 
by (d — w )/(\d — w Q \ + e) 2_p , when p G [1, 2), and letting e — > 0) yields 

(2.39) -4 / Id-tuoPda: 
p dt J M 3 

[u ■ V(d — Wo)](d — u>o)\d — wo\ p ~ 2 dx 



+ I A(d - w )(d - w )\d - w a \ p 2 dx 

2{wq ■{d-w ) + ^\d- w Q \ 2 } 2 \d - w \ p - 2 -hd- w \ p+2 dx 
4 8 

= Ir + h + Ig, 

where the Ig was obtained similarly as in the previous calculation for Iq. We 
estimate If, Is and Ig as follows: 
Integrating by parts over ball Br gives 

h= lim [ / |d — wq\ p u ■ nda — (p — 1) / [u ■ V(d — wo)](d — wo)\d — w$\ p ~ 2 dx\. 

R -^°° JdB R JBn 

It implies that 

(2.40) ph = lim / \d-w \ p u-nda 



R— ► 



9B« 

1/2 / r \ 1/2 



< lim |d-ui | 2 dcr) f / \u\ 2 d(j 

R-S-oo \J d B R J \JdB R 

for any p > 1, where we used that |d — wo| < C. 

By Lemma 12.21 we know that L 3 |d— wo| 2 da; is bounded, and J ffi3 |w| 2 da; is bounded 
from the energy estimate (|2.16|) . Thus, for the inequality ()2.40|) . using arguments 
similar to the ones applied to derive the convergence (|2.27[) , will yield 



/ |d — Wo\ p u ■ nda — > 
J OB R . 



for an appropriate sequence R4 — > 00, for any p > 1. Thus, 

(2.41) Z 7 = 0. 
For Is, integrating by parts over Br yields 

(2.42) 7 8 =lim[/ d(d-w ) ■ (d-w Q )\d-w \ p - 2 da 

R ^°° JdB R on 

\V{d-w Q )\ 2 \d-w a \ p - 2 dx 

B R 

(p-2) [ \\7(d-wo)-(d-wo)\ 2 \d-w \ p - 4 dx} 

J B R 

lim (K1-K2-K3). 

R— >oo 
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The boundary term K\ is estimated as 
(2.43) K\ = 

< I |V(d- w )\\d- u7o| p ~ 1 d<7 



d(d-wp) 2 

■ (d - w )\d - w \' da 

dB R 9n 



dB T , 



< / \\7(d~ w )\da 

JdB R 



where we used that \d — wo\ p _1 < C for any p > 1. 

In Proposition (|2.1j) . let fc = 1, m = 2, r — 1, p = 2 and g = 1. For a = | the 
inequality (|2.18[) yields 

|V(d(*) - w )\dx 

/ r \ 1/7 / f N 5/7 

<C / |A(d(t) — too)|dx J / |d(t)-io |da; 



where we used Ladyzhcnskaya estimate (jl.lOj) and the estimate (jl.llj) . Co is a 
constant depending on initial data and C(t) is the time dependent function in 
(fi~TT|) . Thus, for any fixed t > 0, 



/ |V(d- w )\d<7 -> 



for a sequence i?^ — ► oo. Hence, from (|2.43[) . we have 
(2.44) lim Ki = 0. 

R—tao 

Since 

-|V(d - w )| 2 |d ~ wo\ 2 < -|V(d - iu ) ■ (d - w )| 2 , 

we have 1 

-|V(d - Wo )| 2 M - ^oT 2 < -|V(d - w ) • (d - ™ )| 2 |d - w \ p -\ 
which implies from (|2.42p that 



(2.45) 



-K 2 - K 3 < (p - 1) / -|V(d - w ) • (d - ™ )| 2 |d - w \ p - 4 dx 
Jb r 



4(p-1) 

9 



/ |V|d-wo|^| 2 da;. 

JBr 



Combining p^2]) . (f!Q4"]) and (j2T45|) yields 

(2.46) /8<- 4(P ; 1} / |V|d-™ | 5 | 2 dz 



for any p > 1. Slightly modifying the process to estimate Iq gives 

1 

8rf 



(2.47) h < — I \d-w \ p+2 dx. 
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Combining (f2T39|) with inequalities (|2~4"Tj) . (f2~4"6]) . and (f2~47l) yields 

(2.48) ~ / \d-w Q \ p dx+{p-l) [ \\7(d-w )\ 2 \d-w \ p ~ 2 dx 

P dt J & 3 J R 3 

<C [ \d-w \ p+2 dx. 

JR 3 

Denote v = \d - w \ p/2 - Then (|2"^gj) can be rewritten as 

(2.49) -4- [ \v\ 2 dx+ 4(p ~ 1) / \Vv\ 2 dx<C [ \v\ 2 \d-w \ 2 dx. 

P dt J R 3 v Js.3 J R 3 

Applying Holder inequality and Sobolev inequality to the right hand side yields 

/ r \ 1/3 / f N 2/3 

\v\ 2 \d-w \ 2 dx < [ \v\ 6 dx) / \d-w \ 3 dx 



<C [ \Vv\ 2 dx( [ \d-w \ 2 dx) f \\d~wo\\ 2 L ^ ( 
Jr 3 \Jr 3 / 



From Lemma 12.21 and Lemma 12.31 when p > 1, we can choose initial data small 
enough so that 



(2.50) C^Jd-w \ 2 dx^ ' \\d-w Q fH m 
It follows from (f2~49]) that 

(2.51) ~f \v\ 2 dx+ 2{P ~ 1) [ \Vv\ 2 dx<0. 

P dt J R 3 P Z J R 3 

The inequality (|2 .38[) . for p > 1. can be obtained integrating (|2 . 5 1 1) over [0, T] 
When p = 1, we have from (|2.48|) 

— f | d — wq I dx < C ( \d — wo\ 3 dx 

dt J R 3 J R 3 



< 2(p-l) 



i' 2 



< C / |d - w \ 2 dx\\d - w ||l°=(r3) 

JR 3 

< C [ \Vd\ 2 + \Ad\ 2 dx 

where we used Lemma |2~21 and Lemma |2"31 Integrating over time [0, T] yields 

| d — wo\dx 

<[ \d -w \dx + C [ [ \Vd\ 2 + \Ad\ 2 dxdt < C 
Jr 3 Jo Jr 3 

where we used the energy estimate (|2.16[) . This completes the proof of the lemma. 

□ 



ASYMPTOTIC BEHAVIOR OF SOLUTIONS TO THE LIQUID CRYSTALS SYSTEM IN K 3 13 



2.3. Decay of d — wo and V(rf — Wq). We first establish that d — u>o decays in 
LP(R 3 ), for p > 1 at the rate (1 + t)~^ (1 'p } . 

Theorem 2.6. Let d be the solution obtained in Theorem ll.ll Assume d — wq £ 
L P (R. 3 ), p > 1. Assume the initial data satisfies the conditions of Lemma \2.5\ Then 
for any 1 < p < oo, t > 

(2.52) \\d(;t)-wo\\ LP{m <C(l+tr |(1 ~*\ 

where the constant C depends on X p as defined in Lemma \2.5\ 



Proof: Note that as p — > oo the constants A p in Lemma 12.51 will tend to zero, 
hence we cannot pass to the limit as p — > oo. Therefore this result does not give 
the decay for the L°° norm. 

We proceed by induction for k with p = 2 k . The other powers p follow by interpo- 
lation. When k = the theorem follows by Lemma |2~51 Suppose it holds for s = k, 
then we have 

(2.53) \\d{-,t)-wo\\ I> > 0ia) <C(l + trK 1 -&, 

Let v = \d(-, t) — w \ 2k . 

Recall the inequality (|2.51|) ( which holds provided the data satisfies (|2.50|0 

(2.54) \v\ 2 dx + C p 2(p ~ 1) [ \Vv\ 2 dx < 0, 

dt J H 3 P J & 3 

By Gagliardo-Nircnbcrg we have 

3/5 / . \ 4/5 



J \v\ 2 dx <c(yj \^v\ 2 dx^ (j \v\dx 



Hence using the inductive hypothesis on the last integral on the right hand side we 
have 

- J \\7v\ 2 dx<-C^j \v\ 2 dx^j (l+i) 2 ^- 1 ) 
Combining the last inequality with (|2.54j) yields 

d Jr3 \v\_dx ^ ^ t ->2(2 fc -l) 

*(r R3 H^) 5/3 - 

Integrating and reordering terms yields 

-i 3/2 



\v\ 2 dx < 



.(l + Ci + t)) 2 ^- 1 ^ 1 . 

Since |(2(2 fc - 1) + 1) = |(2 fc+1 - 1) the induction step is obtained, establishing 
the conclusion of the theorem. 

□ 

As a consequence of the last theorem, we derive the decay of V(d — ivq). 

Corollary 2.7. Let d be the solution to system obtained in Theorem 11.11 

Then 

||V(d-tflo)||| a <C(l + t)-t, 
where C depends on initial data. 
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Proof: Take p = 2 in Theorem ESI 
(2.55) \\d(-,t) - w \\ L 2 (R3} < C\\d - w a \\ L 2 {R3) (l + tyi. 

Gagliardo-Nirenberg inequality (|2.18|) yields 

/ \V{d~w )\ 2 dx <C ( [ \d- w \ 2 dx ] ( / | A(d - w a )\ 2 dx ) 

JR 3 \Jr 3 J \JR 3 J 

< C (j \d-w \ 2 dx^j 2 
<C(l + t)-i, 

where in the last two steps we used Ladyzhenskaya energy estimate (|2.17|) and 
(|2.55p . respectively. The constant C depends on initial data. It completes the 
proof. 

□ 



3. Decay of Velocity 

An application of the Fourier Splitting method [53] is used to establish L 2 decay 
of velocity it. 

Theorem 3.1. Let u be the solution obtained in Theorem If additionally 

u Q e L 1 (R 3 ), then 

\\u(;t)\\h<C(l + t)-K 
where C depends on initial data, the L 1 and L 2 norm of Uq. 

Proof: Multiplying the Navier-Stokes equation in system (jl.ll) by u and inte- 
grating by parts yields 

(3.56) / \u\ 2 + 2f \Vu\ 2 dx = 2 [ Vu(Vd®Vd)dx. 

dt JR 3 Jr3 Jr3 

Holder and Cauchy Schwartz inequalities yield 

2/ Vu{Vd®Vd)dx< I \Vu\ 2 dx + C I \Vd®Vd\ 2 dx. 

JR 3 JR 3 JR 3 

Thus, we derive from (|3. 561) 

T / l"| 2 + / Wu\ 2 dx<c( \Vd®Vd\ 2 dx. 

dt Jr3 Jr3 Jr3 

The right hand side of above inequality can be estimated as 

/ \Vd® Vd\ 2 dx = / (Vd®Vd)(Vd(g>Vd)dx 
Jr 3 Jr 3 

= -3 / (d - w ) Ad Vd <g> Vd 
Jr 3 

<\ f |Vd® \7d\ 2 dx + C [ \d-w \ 2 \Ad\ 2 dx, 
2 Jr 3 Jr 3 
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from which it follows 



/ \Vd®Vd\ 2 dx<C [ \d-w \ 2 \Ad\ 2 dx 

<c( f \d-w \ p dx) ( f \Ad\^dx] 
= C\\d-w \\l Pm (^JjAd\ 2+ ^dx 



<C\\d~w \\ 2 LJ 

for p > 2, where the last step followed from Ladyzhenskaya estimate (|2.17p and 
the fact that || Ad|| £ =c( R 3 x [ t jj is bounded since d is regular in the sense stated in 
Theorem ll.il Thus, it follows from Theorem 12.61 that 



|Vd® Vd\ 2 dx < Cil + ty 3 ^-^, 

3 

for any p > 2. Therefore, 

(3.57) ^ / \u\ 2 dx+ f \Vu\ 2 dx<C(l + t)- 3{1 -^. 

dt Jn>3 J R 3 

Applying Plancherel's theorem to (|3.57|) gives 

(3.58) ^/ \u\ 2 d(+ [ |£| 2 H 2 d£<C(l + i)- 3(1 -5). 

dt JR3 J R 3 

The idea is to decompose the frequency domain R 3 in integral J R3 \£,\ 2 \u\ 2 d^ into two 
time-dependent subdomains. The time dependent subdomains are a 3-dimensional 
sphere, S(t), centered at the origin with an appropriate time dependent radius and 
its complement. For this we rewrite p.58|) as 

1 1 \u\ 2 da < - [ \e\u\ 2 dz - 1 \e\u\ 2 dt + c(i + t)-w-$\ 

at JR3 Js(t)" JS(t) 

where S(t) is the ball 

5(t) = {eeM 3 :|e|<r(t) = ( T ^ 7 ) 1 / 2 } 
for a certain fc, which will be determined below. Hence 



1/ \uM<--^-f \&\*dt- 1 \e\u\ 2 dt+c(i+ty 3 ^ 

at ./i»3 i + i J say Js(t) 



'5(t)= JS(t 



N 2 # + / (Az- + c(i + 



1 + 1 V»3 .1 s(t^ 1 + i 



and 



(3.59) ±f \u\ 2 di+-^- ( \uM<-^- f \uM + C(l + t)-^\ 
at Js? 1 + 1 Jm? 1 + 1 Js(t) 

The following estimate, which will be established later, is needed 

(3.60) \&{t,t)\<C\t\- 1 
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for £ € S(t), where C is a constant only depending on the initial data. Combining 
the inequalities (|3.59[) and ()3.60j) yields 



dt I ^ 



l + t 



\u\ 2 d£ < 



C 
l + t 



\Z\- 2 d^ + c(i + t)- 3{1 -^. 



S(t) 



Multiplying by the integrating factor (1 + t) k yields 
d 



dt 



(l+t) k / \ufd£ 



< C(l + 1)^ 1 / iei~ 2 rfe + C(l + t) k - 3{1 ~p\ 



S(t) 



Since p > 2 and 



it follows that 



/ ier a dg < c / : 

JS(t) Jo 



r 2 r- A dr < C(l + t) 



-1/2 



d 

Integrating in time yields 



(l+t) k / \u\ 2 dH 



< C{l + t) k 



(l+t) k \u\ 2 d£< \u(£,0)\ 2 d£ + C[(l + t) 



Thus, 



K 3 



\u\ 2 di < (l + ty k / \u{^0)\ 2 dZ + C[(l + t)-i - (l + t)" 



Since mq € I 2 , it follows that u(0) £ -L 2 from Planchcrcl's theorem. Hence 

H 2 d£<C(l + t)^, 



Hence 



\u\ 2 dx < C(l + t)" 



To complete the proof we need to establish the inequality (|3.60[) . Taking the 
Fourier transform of Navier-Stokes equation in system (jl.ll) yields 

(3.61) u t + \eu = G(Z,t) 
where 

G(f , t) = -T(u ■ Vu) - J"(Vp) - J"(V • (Vd <g> Vd)), 
and J 7 indicates the Fourier transform. Multiplying (|3.61[) by the integrating factor 
e^l 4 yields 

Integrating in time gives 

(3.62) u(£,t) = e- | « |2 *w + / e~ l « |2( *- s) G(e,s)ds. 



We assume for the moment the following auxiliary estimate, which we will prove 
below, 



(3.63) 



\G(U)\ <C\£\ 
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Combining (f3~62]) and (fX63|) yields 

(3.64) \u(£,t)\ < e-'«' 2t KI + f e-^ 2{t - s] \t;\d S . 

Jo 

Since uq gL 1 , we have |uo| < C for all £ and some constant C. Performing integra- 
tion in (|3.64[) gives 

m,t)\ < ce-i^i 2 * + ^(i - e -i«i 2 *) < cier 1 

for £ G S'(i). To finish the proof we need to establish (|3.63|) . For this purpose we 
analyze each term in G(£,t) separately. We have 

\T{u- Vu)| = |.F(V- (u®u))\ < V / |uV||<j|dar. 

Since u G L°°(L 2 ) by the basic energy estimate, we have 

\F{u- Vu)\ < C|£|. 

By the basic energy inequalities (|2.161 I2.17[) we have Vd G L°°(L 2 ) proceeding 
similarly as for the last inequality we have 

|.F(V-(Vd®Vd))| < C|£|. 
Taking divergence of Navier- Stokes equation in system gives that 

Taking the Fourier transform then yields 

\t\ 2 F(p) = -£)&&.F(uV) - ^^-^(vrf'v^). 

Since .F(itV) G L°° and J"(Vef V<#) G it follows that 

J-(p) < C, 

and thus J-(\7p) < C\£\. It completes the proof of (|3.63p and hence completes the 
proof of theorem. 

□ 

Combining Theorem 12.61 Corollary 12.71 and Theorem 13.11 yields the proof of the 
main Theorem 11.21 

□ 

Remark 3.2. The decay rate for the velocity u in L 2 obtained in [35], for the 
bounded domain case, is (1 + t)~ 1 - 2 » where 6 G (0, 4). When 9 is close to 0, 
then — 1 ^ 2 g would be very small, meaning the decay is very slow. In this paper, 

we obtained the decay rate for velocity u in L 2 with (1 + a fixed constant 

algebraic rate. The advantage comes from the fact that we work on the whole space 
R 3 where we can apply the Fourier splitting method. 
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Remark 3.3. It was pointed out in the first section that there is an essential 
difficulty to apply Lojasiewicz-Simon approach in whole space R 3 . However, in 
weighted Sobolcv spaces of M 3 , the compactness is recovered. Thus, wc expect 
there is hope to construct certain Lojasicwicz-Simon type inequality in weighted 
Sobolev spaces and proceed with the method in [28] to derive the decay of solutions 
to the LCD system in weighted Sobolcv spaces. 



Appendix A. Existence of Classical Solutions in M 3 

In this section we sketch a brief proof of the existence Theorem 11,11 Section 1 . 
As mentioned in the introduction, for bounded domains in R 3 , the existence of 
global regular solutions to the flow of nematic liquid crystals with constant density 
has been established in [19] provided the viscosity is large enough. The existence 
of global regular solutions to the flow of nematic liquid crystals with non-constant 
density has been established in [6] provided the initial data is small enough. In 
both of the above papers, a Ladyzhenskaya energy estimate (higher order deriva- 
tive estimate) was derived and hence a relatively standard bootstrapping argument 
yielded a regular solution. 

The proof of Theorem 1 1 . 1 1 will be given through four steps. In the hrst step, on 
a sequence of balls B Rn with radius R n , centered at the origin, we obtain the ex- 
istence of a Galcrkin approximated solution (u n ' m ,d n,m ) for the system (|1.1[) with 
modified initial data, for each m = 1,2, 3, .... In the second step, we establish an 
estimate of d n,m — wo in L 1 (_B/f ri ) for any fixed time t > 0. In the third step, we 
take the limit m — > oo. In the forth step, we take the limit R n — > oo. In fact, we 
are able to show that all the estimates in Ladyzhenskaya energy method in step one 
are independent of the domain size. Thus we can take a subsequence of solutions 
on balls B Rri which converge to a limit in M 3 when R n goes to infinity. 

Lemma A. 1. Assume u a e H 1 (R 3 ) and d = d (x) - w E H 2 (R 3 ) n L 1 (R 3 ) with 
\wq\ = 1. There exists a sequence of functions {(ug , oIq )}^° =1 and a sequence of real 
numbers {Rn\^ = i with R n — > oo as n — > oo such that, 

(A.65) v% e Hl(B Rn ), with itg — > uq in H (M ) as n — > oo 

(A.66) J l € H$(B Rn ), with <% -» d in H 2 (M. 3 ) as n -> oo 

where d^ = c?g — wq . Moreover, 

(A.67) \K\\ LH B Rn} < C\\d Q \\ LHm , 

and \d%\ < 1. 

Proof: Such a sequence of functions can be constructed easily as follows. Let 
be a sequence of smooth functions such that 

Cn(a:) = \ 

[0, m B 2 R n 
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and \Cn(x)\ < 1 for all x G R 3 . 

Define Uq = Cn^o an d dg = ( n do. Let d^ = d' a l + wq. These £„ can be chosen such 
that (j A.65|) . (|A.66[) and (|A.67[) are satisfied. In addition, we have 

Kl = \(nd + (i - Cn)H < Cn + (1 - (n) = 1- 

□ 

For the sequel we assume as initial conditions 
( u(x, 0) = Uq(x), 

(A ' 68) {,i(x,o) = «M,' i ° S ' i " x[0 ' n 

where Ug and dp are as obtained in Lemma I A. 1| and the boundary conditions 

f u(x,t)=0, 
(A.69) on9i? fl „ x [0,T). 

I d(x, t) — wo = (J, 

We have the following existence of solutions to system with (|A.68j) and 

(|A.69j) on each ball B Rn . 

Theorem A. 2. Let B Rn be the ball centered at origin with radius R n in R 3 . As- 
sume that uq G ii^R 3 ) and do — wq G H 2 (M. 3 ). The system U.l\) with initial and 
boundary conditions \A . 68\) and HA . 69\) has a smooth solution (u n ' m ,p n ' m ,d n ' m ) 
for each m = 1, 2, 3, ... satisfying, for any T > 

u n,m g L 2 (0, T; Hi (B R J) n L°° (0, T; L 2 (B Rn )) 

JH.m g L2( 0j r;ir2( BflB )) n£ ao( 0jr;ir l( BjjB jj 

with d n,m = rf™' m — hjq and |d n,m | < 1. Moreover, it satisfies the energy inequality 
(A.70) ^ / i|u"' m | 2 + i|Vd"' m | 2 + F(d n < m )dx 

+ / |Vu"^ l | 2 + |Ad"^ l -/(d n ' m )| 2 dx<0. 

Proof: The existence proof is obtained through the standard Galerkin approx- 
imation method, see [115] and [5] ■ We only need to give a brief explanation on the 
claim d n,m | < 1 by applying a maximum principle argument. Notice that the ap- 
proximated initial data dp satisfies |dg(:r)| < 1 for all i£l 3 . Suppose there exists 
a point (xo,to) in the interior of the domain B Rn x [0,T), such that |d n ' m | 2 attains 
a maximum value at this point. Multiplying the equation 

d n, m + u n, m _ vd n., m = Ad n,m _ J_ (| d n,m|2 _ ^n.m 

n 2 

by d n ' m yields 

(A.71) ^l d " : '"| 2 + u "' m ' V l d "' m | 2 ' dn '" 1 



A|d"' m | 2 - 2|Vd n - m | 2 - -A-(|d"< m | 2 - l)|d"< m | 2 . 



2 

rf 

At the maximum point (x ,t ), we have £ t \d n < m \ 2 = V|d"> m | 2 = and A|d"' m | 2 < 
0. Thus, it follows from the equation (|A.71[) that, at the point (xo,to) 

(\d n ' m \ 2 - l)|d"< m | 2 < 0. 
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This insures that |d n,m | < 1 at any interior maximum point (xo,to). Therefore, 
\d n > m \ < 1 for all points in B Rn x [0, T). 

□ 

For the solution (u n,m ,d n,m ) obtained in the above theorem on B Rn x [0,T), we 
define energy quantity 

< m (t) = ||V« W ' ro ]|i 2(SRj + \\A(d^ m Wo)\\l HBRn y 

With a slight modification of the proof of Theorem 3.1 in [6], we are able to show 
that 

Theorem A. 3. Assume that uq G H 1 (Bn ri ) and da G H 2 (Bn n ), and ||«o||#i(R3) + 
||<io|| jj2flg3) < 00 ■ Let (u n > m ,d n ' m ) be solutions obtained in Theorem \A.2\ There is 
a positive small number eo such that if 

(A.72) IM|hi( R 3) + ||do - wo||^ 2 ( R3) < e , 

then 

(A.73) f \Vu n - m \ 2 + \A{d n - m ~w )\ 2 dx 

+ f [ \Au n - m \ 2 + \VA(d n ' m -w )\ 2 dx 

Jo J B Rn 

< G{\\u \\ 2 H i {BRn) + \\d - w \\ 2 H 2 {BRn) ), 

for any T > 0, where the constant C is independent of domain size R n and m. 

There is no need to prove the theorem except that we need a brief explanation 
on the last claim that constant C is independent of R n . In the proof of Ladyzhcn- 
skaya energy estimate in [6], we only use the Gagliardo-Nirenberg interpolation 
inequalities and standard elliptic inequalities. That is we use 

K'lll* <C\\u n ' m \\ L 4\7u n - m \\l 2 

\\\7(d n - m - w )\\U < C\\V(d n - m - w )\\ L 2\\A(d n - m - w )\\h 
and the elliptic estimate 

\\D 2 u n > m \\ L 2 < C\\Au n - m \\ L 2 

for u n - m and d n ' m — wo vanishing on the boundary. In the above inequalities, the 
various constants C are independent of the size of the domain. Thus the constant 
C in (|A.73|) is independent of R n . 

For the second step, we derive a time dependent estimate of d n,m (t) — wq in 
L\B Rn ). 

Lemma A. 4. Let d n,m be the solution obtained in Theorem \A.HA In addition, 
assume do — wq G L 1 (M 3 ). Then 

(A.74) f \d n > m {t)-w \dx<{C t+ [ \d - w \dx)e ct 

JB Rn JR 3 

where the constant C'o only depends on initial data and constant C only depends on 
V- 
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Proof: By the second equation in we have, 



(A.75) 
d_ 
dt 



4/ \d n < m (t)-w \dx 



III .III 



(d n > m (t)-w )-d? ^ 

\d n > m {t)-w \ X 

(d n - m (t) - w ) ■ Ad n > m (t) , f (d n - m (t) - w ) ■ [u n - m ■ Vd n > m (t)] , 

dx — I r~. ; dx 



BRn \d n >™(t)-w \ J BRn \dP"™{t)-W \ 

(d n < m (t)-w a )- f(d n > m (t)) ^ 
\d^(t)-w \ 



h+h+ I: 



There is no need to worry about the singular points of (d n ' m — wq)^ 1 in the above 
equation, since each term on the right hand side contains ^„' m ~^°| . We deal with 
the three terms I\ , I2 and I3 in the following way. 

Since d n,m (t) — wo = on the boundary dB^ n . Replacing dt by its value in the 
direction equation and integration by parts yields 

(a.76) h =-[ wfy-^\ 2 dx 

+ f \(d n ' m (t) - w ) ■ V{d n - m {t) - w )\ 2 \d n > m {t) - wo\- 3 dx 

JB Rn 

<0, 

since \(d n > m (t) - w ) ■ V{d n ' m {t) - w )\ 2 < \{d n ' m (t) - w )\ 2 \V{d n ' m (t) - u> )| 2 . 
By Ho ldcr inequality we have 

(A.77) \h\< [ \u n > m \\V(d n > m -wo)\dx 

JBRn 

\u nm \ 2 dx\ ( I \Vd n ' m \ 2 dx\ <Cn, 



where we used the energy estimate (|A.70[) . and the constant Co only depends on 
the initial data. 

Recall that by definition f{d n > m ) = ^{\d n ' m \ 2 - l)d"< m , and |d"' m | < 1 from 
Theorem I A. 2 1 and |wo| = 1, hence 

(A.78) \I a \ < 1 / \f(d n ' m )\dx 

<C I \d n ' m - wo\\d n ' m + w \\d n ' m \dx <C [ \d n ' m - w \dx, 

where the constant C depends on r\. Combining the inequalities (|A.75[) . (|A.76I) . 
(|A~77l) and (|A~78l) yields 



^ I \d n - m (t) -w \dx < C I \d n < m {t) -w \dx + C . 

Bn n 



(It 
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Integrating over [0,t], and Gronwall's inequality (see [10]) gives, 

f \d n - m (t) - w \dx < (C t + J |dff - w \dx)e ct , 

for any t > 0. Inequality (|A.74[) now follows from the last estimate in Lemma I A. II 
This concludes the second step. 

In the third step, we take the limit of the Galerkin approximating solutions 
(w"' m ,cT< m ) as to -> co. By the estimates (jA~70| and (|A"7T3jl . there exists (u n ,d n ) 
for each n = 1, 2, 3, ... such that, taking subsequence if necessary, 

u n..m _± u n weakly m jr. #1(£ R J) ; 

u n,m u n strongly m L°° {{) ,T ' L 2 [B R J) , 

d n, m ^ d n W eakly in L 2 (0,T;H 2 (S H J), 

d n,m _^ d n strongly m jr. H ^ Br J) with rf ™ = Wq on 

It follows easily from the above convergence that (u n ,d n ) is a weak solution to 
the system (|1.1[) with initial condition (|A.68[) and boundary condition (|A.69[) on 
Br„ x [0,T). Moreover, the solutions (u n ,d n ) satisfy the basic energy inequality 

(A.79) l K f + l ]Vd nf +F{d n )dx 

+ [ \\7u n \ 2 + \Ad n - f(d n )\ 2 dx <0, 
and the higher order energy inequality 
(A.80) / |Vu"| 2 + |A(d"-u;o)| 2 da; 

+ / f \Au n \ 2 + \VA(d n -w )\ 2 dx 

JO JB Bn 

< C(\\Uo\\ 2 H i(B R J + No - w o\\H2 {BRn) ), 

for any t > 0, where the constant C is independent of domain size R n . In addition, 
from Lemma IA.4I it follows that d n satisfies the estimate 

(A.81) / \d n {t)-w \dx<(C t+ [ \d - w \dx)e ct 

JB Rn JK 3 

where the constant Cq only depends on initial data and constant C only depends 
on rj. 

In the forth step, we extend the solutions (u n ,d n ) on B Rn to the whole space 
R 3 by taking limit R n — > oo. With the estimates (|A.79|) and (|A.80|) we can extract 
a subsequence {(u lk , d lk )} k from (u n ,d n ) for n > 1 such that 

u lk ^u {1) in L 2 (0,T;H\B Rl )) 

u lk ^u (1) inL°°(0,T;L 2 (B Kl )) 

d lk ^d (1) inL 2 {0,T-H 2 {B Rl )) 

d lk ->d (1) in L 2 {Q 1 T ] H l (B Rl )) 
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and the limit (um, d(i)) satisfies system (jl.ip in distribution sense and the estimates 
(TP!) (1X80] on B Rl x [0,T). 

On _Br 2 x [0,T), we take subsequence {(u 2fe , d 2k )}™ =1 from {(u lfc , ^ lfc )}^ 1 such 
that {u 2k }°^_ v {p 2k }^_ 1 and {d 2 * 1 }?^ converge to U( 2 ) and d(2) respectively in 
the same convergence sense as above. And we have that 

u (2)\b x = W(l), ^(2)|-Bi = ^(l)- 

Repeating the process on each x [0, T), we can take subsequence {(u nk , c£ nfc )}^ 1 
from the sequence {(u(™" 1)fc , d ( "~ 1)fe )}^ 1 , such that {" nfe }^ 1 and {d™' }^ con- 
verge to it(„) and d(„) respectively. And we have that 

"(n)|s„_! = 1t( n -l), d(n)|s„_i =d(„-l). 

Then we take the diagonal sequence {(u kk , d kk )} € ^'_ 1 and let k —> oo. This sequence 
(if necessary, take a subsequence of it) converges to (u, d), in M 3 x [0, T). The limit 
(u, c?) satisfies the system fll.l[) in the sense of distributions and satisfies the energy 
estimates 

(A.82) f \u\ 2 + \Vd\ 2 + 2F(d)dx + 2 [ [ \Vu\ 2 + \Ad - f{d)\ 2 dxdt 
Jr 3 Jo Jr 3 

< [ \u \ 2 + \Vd \ 2 dx 



(A.83) 



/ \Vu\ 2 + \Ad\ 2 dx + f [ \Au\ 2 + \\7Ad\ 2 dxdt 
Jr 3 Jo Jr 3 



< c(I|woIIhi(r3) + IMo - woIIh2 (r3) ). 

In addition, the solution d satisfies the estimate 

\d{t) - w \dx < (C t + | do -w \dx)e Ct 
Jr 3 

with constant C'o only depending on initial data and constant C only depending on 
V- 

Estimates (|A.82[) and (|A.83[) allow us to apply the "bootstrapping argument" 
as used in [6] and [19] , and prove that the limit (u, p, d) is a classical solution to 
system satisfying the desired estimates in Theorem 11.11 This completes the 
proof of Theorem 11.11 
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